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Abstract. We show how a large class of sufficient conditions for the existence 
of bound states, in non-positive central potentials, can be constructed. These 
sufficient conditions yield upper limits on the critical value, g^ , of the coupling 
constant (strength), g, of the potential, V(r) = —gv(r), for which a first ^-wave 
bound state appears. These upper limits are significantly more stringent than 
hitherto known results. 
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1. Introduction 

There exist in the literature several necessary conditions for the existence of at least 
one £-wave bound state in a given central potential. These necessary conditions yield 
lower limits on the critical value, , of the coupling constant (strength), g, of the 
potential, V(r) = —gv(r), for which a first ^-wave bound state appears. 

In 1976, Glaser et al. have obtained a strong necessary condition for the existence 
of bound states in an arbitrary central potential in three dimensions (h 2 /(2m) = 1) 

HI 

(p~lf^T(2p) [°°dr l2 F 

where V~(r) = max(0, — V(r)) is the negative part of the potential and with the 
restriction p > 1. This inequality is nontrivial provided that the potential V(r) is less 
singular than the inverse square radius at the origin and that it vanishes asymptotically 
faster than the inverse square radius, say (for some positive e) 

lim [r 2 ~ e V(r)] =0, (2) 
lim \r 2+£ V(r)} = 0. (3) 

r — >oo 

We assume throughout that the potentials satisfy the relations © and © and that 
they are piecewise continuous for r € ]0,oo[. The lower limit on g^ obtained from 
is actually very accurate as it has been demonstrated on several examples (see for 
example PEIE] as well as Section |3}. 

Recently other strong necessary conditions have also been obtained [3] 

2 

(27+Tp 

6 

(2* + 1)3 

As shown in [5J, these two inequalities, J3J and J5J, are natural extensions of the 
Bargmann-Schwinger necessary condition 0]O3 (first obtained by Jost and Pais [H]) 

1 f°° 

dxxV~(x)>l. (6) 



dxx~ M V-(x) I dyy zt+z V-(y)>l, (4) 
dxx- 2e V-(x) I dyyV-(y) j V dzz 2e+2 V-{z)>l. (5) 



21+ 1 

Actually the inequalities 10, Q and (J5J are the first members of a sequence of 
necessary conditions which yield a monotonic sequence of lower limits on the critical 
value of the strength of the potential, gi^, which converges to the exact critical 
strength [3]. This remark implies that the inequality (0 yields stronger restriction 
than the relation @. The complexity of each member of this sequence of necessary 
conditions becomes rapidly important and only the relation 10} and (JSJ can be easily 
used. It has been shown, with some test potentials, that the relation JSJ can be better 
than the relation Q, especially for i — (see tests performed in Ref. |3] and in 
Section below). 

Other necessary conditions for the existence of bound states can be found in the 
literature (see for example 0IH] and for reviews see |9l ll()llll| ). but none, in general, 
yields stronger restrictions than (JJJ and (j^J. 
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Few sufficient conditions for the existence of a ^-wave bound state in a central 
potential, yielding upper limits on gc , can be found in the literature. Let us mention 
two sufficient conditions found by Calogero in 1965 ^1 ITS'] 

/>a />oo 

/ drr\V(r)\ (r/af e+1 + / dr r \V{r)\ (r/a)^ +1 ) > 2£ + 1, (7) 

JO J a 

and 

r°° -i 

a dr\V(r)\[(r/a) 2e + {r/a)- 2e a 2 \V(r)\] >1. (8) 
Jo 

These two conditions apply provided the potential is nowhere positive, V(r) = 
— |V(r)|; in both of them a is an arbitrary positive constant, and of course the most 
stringent conditions obtain by minimizing the left-hand sides of Q and (jSJ) over all 
positive values of a. 

Few other sufficient conditions for the existence of bound states can be found in 
the literature (see [21 EH El); but they are either quite complicated or less stringent 
than and ©. 

In this article, we obtain a strong sufficient condition for the existence of 
bound states yielding accurate restrictions on the critical strength g^ which improve 
significantly the restrictions provided by the relations (7J and ©■ 



2. Sufficient condition and upper limit on the critical strength 

(£) 

The idea used to derive the upper limit on gc is to transform the standard 
eigenvalue problem obtained with the time independent Schrodinger equation, where 
the eigenvalues are the eigenenergies, into an eigenvalue problem where the eigenvalues 
are the critical coupling constants. These critical values of the strength of the potential 
correspond to the occurrence of an eigenstate with a vanishing energy. 

Following Schwinger jS] (see also |15p. we consider the zero energy Schrodinger 
equation that we write into the form of an integral equation incorporating the 
boundary conditions 

/•OO 

u e {r) = - dr' g e {r, r') V{r') u t {r'), (9) 



JO 

where ge(r, r') is the Green's function of the kinetic energy operator and is explicitely 
given by 

^ r > r ') = ^TT r < +lr >"' (10) 

where r< = min[r, r'\ and r> = max[r, r']. An important technical difficulty appears if 
the potential possesses some changes of sign (see relation (|1 1|) below) . This is overcome 
in the derivation of necessary conditions, or of upper bound on the number of bound 
states, by considering the negative part of the potential instead of the potential itself 
(V(r) — > V~(r) — max(0, — V(r))). Indeed, the potential V~{r) is more negative than 
V(r) and thus a necessary condition for existence of a £-wave bound state in V~(r) 
is certainly a valid necessary condition for V(r). This procedure can no longer be 
used to obtain sufficient conditions. For this reason we consider potentials that are 
nowhere positive, V(r) = —gv(r), with v(r) > 0. To obtain a symmetrical kernel we 
now introduce a new wave function as 

Ur) = \V{r)\ 1 '*u l {r). (11) 
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Mr) =9 / dr' K e {r,r')Mr'), (12) 
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where the symmetric kernel Ke(r, r') is given by 

K £ (r, r') = vir) 1 ' 2 g t (r, r>) v(r')^ 2 . (13) 

The relation (|12fl is thus an eigenvalue problem and, for each value of £, the smallest 

(a) 

characteristic number is just the critical value g^ . The other characteristic numbers 
correspond to the critical values of the strength for which a second, a third, £-wave 
bound state appears. The kernel (|13l) acting on the Hilbert space L 2 (M.) is an Hilbert- 
Schmidt operator for the class of potentials defined by and (J3J). Thus this kernel 
satisfies the inequality 

'■OO 

dxdy K e (x,y)K e (x 1 y) < oo. (14) 

Consequently the eigenvalue problem (|12(1 always possesses at least one characteristic 
number |16l pp. 102-106] (in general, this problem has an infinity of characteristic 
numbers). Note also that the kernel l|13fl is the so-called Birman-Schwinger kernel 

0E]. 

Now we use the theorem (see for example PP- 118-119]) which states that, 
for a symmetric (positive) Hilbert-Schmidt kernel, we have the variational principle 

max[/ / dxdyKi(x,y)tp(x)ip(y)] = - 1 pr, (15) 
f Jo Jo g ( c> 

for ip(r) satisfying 

poo 

drip(r) 2 = l. (16) 



n 



The maximal value is reached for ip(x) — <f>1{x), where 4>2(x) is the eigenfunction 
associated to gc ■ Consequently for an arbitrary normalized function, f(x), we obtain 
the following upper limit on g c ^ 



poo poo 



dxdyK e (x,y)f(x)f(y) . (17) 

10 Jo 

To apply the above theorem, we simply choose 

f(r) = A[r 2p - 1 v(r) p ] 1/2 , p > 0, (18) 
where A is a normalization factor. With the choice (|18|) . the upper limit (|17[) reads 



(19) 



p oo r poo px 

#<£/ dxF(2p-l;x) / dxF(p;x)x- c dyF(p;y)y L 
Jo Uo Jo 

with F(q; x) = x q v{x) i - q+1 ^ 2 and L = £ + 1/2. 

We do not consider other choices for the function f(r) here since, as shown in 
Section [21 the relation (|19|l is already very accurate. We just mention that another 
possible choice for monotonic potentials is f(r) — A[v(r)(v(0) — v(r)) p } 1 ^ 2 . We have 
verified with an exponential potential, see 11221) below, that this choice yields a slight 
improvement. 
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Table 1. Comparison between the exact values of the critical coupling constant 
<7c ' of the square well potential 1211 for various values of £ and the lower limits 
on obtained with the relations Q, J^J and J§J, called respectively 9qq mt , 
<7g' and g^g and the upper limits obtained with the formula, J7}, JSJ and 1191 . 
called respectively, g£\, g£\ and g{^ w . 



£ 


W 
9bs 


(A) 
9 b 


09 
5ggmt 


yc 


09 

3Ncw 


09 
9ci 


09 

9 C 2 





2 


2.4662 


2.3593 


2.4674 


2.4747 


2.6667 


4 


1 


6 


9.8132 


9.1220 


9.8696 


9.9934 


11.719 


10.068 


2 


10 


19.895 


18.454 


20.191 


20.604 


25.413 


20.895 


3 


14 


32.383 


30.245 


33.217 


34.099 


43.570 


35.424 


4 


18 


47.064 


44.425 


48.831 


50.357 


66.089 


53.519 


■5 


22 


63.788 


60.947 


66.954 


69.295 


92.909 


75.114 



Table 2. Same as for Table HI but for the exponential potential 1221 . In the 
column p, we report the values of the variational parameter p which optimize the 
upper limit fl9l . 



1 


09 
Sbs 


09 


09 
#GGMT 


a {l) 
yc 


SNew 


W 
9ci 


(£) 

9 C 2 


P 





1 


1.4422 


1.4383 


1.4458 


1.4467 


1.6755 


1.5442 


1.4686 


1 


3 


6.8546 


7.0232 


7.0491 


7.0584 


9.7188 


7.7262 


2.4313 


2 


5 


15.257 


16.277 


16.313 


16.334 


24.724 


19.794 


3.4103 


3 


7 


26.265 


29.218 


29.259 


29.289 


46.985 


37.791 


4.4015 


4 


9 


39.616 


45.849 


45.893 


45.932 


76.586 


61.758 


5.3874 


5 


11 


55.120 


66.173 


66.219 


66.264 


113.55 


91.708 


6.3804 



The upper limit ljl9|l is optimal in the sense that it can be saturated, for p = 1, 
by a Dirac-delta potential, V(r) = —g5(r — R), which admits a bound state as soon 
as g = R" 1 . 

Obviously, the sufficient condition for the existence of a f-wave bound state, from 
which the upper limit (|19fl on is obtained, reads 

I* OO f-X f poo ^ — 1 

dxF{p-x)x- c dyF{p-y)y c \c J dxF(2p - 1; x)\ > 1, (20) 
with F(q; x) = x q |V"(a;)|(« +1 )/ 2 , C = £ + 1/2 and p > 0. 
3. Tests 

In this Section, we propose to test the accuracy of the upper limit l|19|) with four 
potentials: a square well potential, 

V(r) = -gR- 2 6(l-r/R); (21) 

an exponential potential 

V(r) = -gR- 2 exp(-r/i?); (22) 
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Table 3. Same as for Table HI but for the Yukawa potential 1231 . In the column 
p, we report the values of the variational parameter p which optimize the upper 

limit rnh . 



£ 


W 


W 
9 b 


5ggmt 


i/c 


5Ncw 


9ci 


(£) 

9 C 2 


P 





1 


1.6689 


1.6643 


1.6798 


1.6826 


2.0505 


1.6810 


1.7217 


1 


3 


8.5999 


9.0384 


9.0820 


9.1039 


13.390 


10.706 


3.1281 


2 


5 


19.553 


21.839 


21.895 


21.937 


35.255 


28.374 


4.5302 


3 


7 


33.931 


40.074 


40.136 


40.194 


67.914 


54.819 


5.9344 


4 


9 


51.368 


63.744 


63.809 


63.880 


111.42 


90.071 


7.3404 


5 


11 


71.615 


92.850 


92.918 


92.998 


165.80 


134.14 


8.7481 



a Yukawa potential 

V(r) = -girR)- 1 exp(-r/R); (23) 
and the STIS (Shifted Truncated Inverse Square) potential 

V(r) = -g(R + ry 2 for < r < aR 

= for r > aR. (24) 

In these potentials, the radius R is arbitrary (but positive) and a is an arbitrary 
positive number. 

The minimization of the upper limit i|19|l over the positive values of p can be 
performed analytically only for the square well potential. We find 

# <c(VC + T+l) 2 . (25) 



The comparisons between the exact value of the critical coupling constants of the 

(£) . .... rn 

potentials, gc , the previously known upper and lower limits reported in Section [T] 
and the new upper limit (|19fl is given in the Tables [21 and [31 for various values of I 
and for the potentials H21fl - H23|) . These comparisons show clearly that the new upper 
limit is very cogent as well as the lower limit {Q obtained by Glaser et al. We have 
also performed other tests, that we do not report here, with nonmonotonic potentials 
and the results obtained are quite similar to those reported in these Tables. 

In Table 01 we present the same comparison for the STIS potential but for I = 0. 
For this potential, the critical coupling constant depends on a. The value of gi°^ is 
obtained, for a given a, by solving the following equation |10| 

A ln(l + a) + 2 arctan(A) = 2tt, (26) 



with A = \J 4 y , c°' ) — 1- For all values of a the results obtained with the new upper 
limit are again very stringent compared to previously known limits. 

4. Conclusions 

The sufficient condition l|20|) proposed in this article yields the upper limit (|19fl on gi^ 
which is analogous to the lower limit obtained three decades ago by Glaser et al. P^. 
The upper limit apply provided that the potential is nowhere positive and that it is less 
singular than the inverse square radius at the origin and that it vanishes asymptotically 
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Table 4. Same as for Table |I] but for the STIS potential J^J and 1 = 0. In the 

column p, we report the values of the variational parameter p which optimize the 
upper limit (TH). 



a 


(0) 

Sbs 


(0) 

9 b 


(0) 

Sggmt 


o (0) 
9c 


(0) 
5Ncw 


o (0) 

yci 


o (0) 

9C2 


P 


0.1 


227.22 


282.11 


269.84 


282.26 


283.12 


306.01 


440.67 


1.2329 


0.5 


13.864 


17.613 


16.842 


17.626 


17.683 


19.311 


24.664 


1.2608 


1 


5.1774 


6.7253 


6.4307 


6.7319 


6.7550 


7.4520 


8.6588 


1.2889 


5 


1.0434 


1.4837 


1.4214 


1.4875 


1.4939 


1.7201 


1.5799 


1.4159 


10 


0.67168 


1.0066 


0.96638 


1.0107 


1.0156 


1.1998 


1.0304 


1.5004 


50 


0.33882 


0.58085 


0.56233 


0.58684 


0.59085 


0.74673 


0.59855 


1.7633 



faster than the inverse square radius. We could use the method proposed in Ref. [T7] 
to consider potentials with some positive parts but the result would then be much less 
neat and then less interesting. 

The method we use to derive the upper limit on the critical strength g^ is 
quite general and other (possibly more complicated) families of upper limits yielding 
(possibly) stronger restrictions on g c ^ could also be obtained. Indeed, the method is 
based on a variational principle for which a trial zero energy wave function is needed. 
There is no limitation on the accuracy of such a trial function, which imply that there 
is, in principle, no limitation on the accuracy of the upper limit on derived with 
this procedure. In this article we have proposed in Section [2] a compromise between 
accuracy and simplicity of the final formula. The accuracy of the upper limit on g^ 
was then tested in SectionElwith some typical potentials. Clearly, the upper limit (fTTjf) 
proposed in this article improves significantly the restriction on the possible values of 
gc obtained with previously known upper limits. 
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